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ABSTRACT: This research work aims to analyse isotropic triangular mild steel plate constrained by boundary 
condition using Finite Element Analysis with CST as master element. Static deformation analysis of a mild steel 
triangular plate fixed at one end under point loads at the free end has been performed using Constant Strain Triangle ( CST) 
as finite element. FEA analysis techniques are used to find the deformation of the plate. Problem Modelling and simulation 
of the plate is done in ANSYS and the results obtained by the finite element analysis using Constant Strain Triangle ( CST) 
are compared. Good agreement has been found out in the deformation results by both the methods. Simulation results are 
critically studied and salient conclusions have been drawn. 
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I. INTRODUCTION 

Finite element analysis (FEA) is a powerful computational technique used for solving engineering problems having 
complex geometries that are subjected to general boundary conditions. In the FE analysis technique, the system is divided 
into discrete finite number of parts or elements by expressing the unknown field variable in terms of the interpolation 
functions within each element. The ease of computation depends on the linearity and the nature of the interpolation 
functions. This paper focuses on the use of the three node triangular element as the finite element. Geometry and 
displacement variables follow isoparametric representation . It is also known as Constant strain triangle(CST) because the 
localized strain in the element remains constant as it depends only on the nodal displacements. The shape functions for this 
elements allow ease of calculations. The paper first relates the shape functions to the displacement variables and element 
stiffness matrix for the CST element has been found out. Later the deformation in the plate is studied by finding out the 
nodal displacements in given loading and boundary condition. Accuracy of the results have been found out by comparing the 
results with ANSYS analysis. 

II. DESCRIPTION 

A triangular mild steel plate is considered in this paper for analysis. Mild steel is the most common form of steel 
used in industries because of its low cost and good mechanical properties. Mild steel is a general term for a range of low 
carbon steel ( about 0.3%). They have good strength and can be bent, worked or can be welded into an variety of shapes for 
uses from vehicles to building materials. A typical CAD model for such a plate is shown in FIG .1 The model is generated in 
the designing software package PRO-E. 




FIG. 1 



For the deformation analysis a typical boundary condition is considered. One edge of the plate is fixed and point 
loads applied at the free end of the plate. The objective is to study the deformation in the plate under the influence of the 
loads that is the amount of deflection of the free end both in horizontal and vertical direction. For the analysis Finite element 
analysis method is utilized and Constant strain triangle is taken as the basic element type. For the purpose of analysis typical 
values of the point loads are assumed and the condition of plane stress condition is considered. 
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2.1) Constant Strain Triangle (CST) : 

In the finite element method, the displacements at points inside the element need to be represented in terms of the 
nodal displacements of the elements. For the constant strain triangle, the shape functions or the interpolation functions are 
linear over the element. The independent shape functions are represented by and r\ . CST will have three nodes at each of 
the vertices of the triangle and each having two degrees of freedom. Let us define the DOF in X positive X axis as q 2il and 

DOF in Y axis as q 2i . where, I is the node number. FIG. 2 shows the three nodes 1, 2 and 3. Thus q x , q 2 , q 3 , q 4 , q 5 , 

q 6 represent the displacements at the three nodes. 
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FIG. 2 
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The displacements inside the element are now written using the shape functions and the nodal values of the 
unknown displacement field. 

u = (4i-4 5 )^+(q3-ft>7+45 

v = (^ 2 -^ 6 )^ + (q4-%) 7 7 + ^6 

For the triangular element, the coordinates x, y can also be represented in terms of nodal coordinates using the same 
shape functions. This is isoparametric representation. From reference [1], for the triangular element, 

C = Bq (1) 

Where, ^ is the strain matrix and q is the displacement matrix. B is a (3 x 6) element strain - displacement matrix relating 

the three strains to the six nodal displacements and is given by 



B = 0 x 32 0 x 13 0 x 21 (3) 

det 7 

Where, X tj 
i and j represent the node numbers. 

det J is the Jacobian of the transformation. Analysis shows that Value of det J is equal to twice the area of the triangular 
element [2]. 

det J = 2 x A e (4) 

Where A e is the area of the triangular element. 

Thus it can be seen from equation (1) that the strain in this triangular element remains constant throughout and that 
is why a lot of computational time is saved by using Constant strain triangle as model element when using Finite element 
method. This paper aims at verifying the computational accuracy of this method when compared to the ANSYS analysis 
method. 

The element stiffness matrix k e is given by 
k e = tAB T DB 

e e 

Where, 

t e = Thickness of the element 

A e = Area of the triangular element 

D = Material property matrix 

Value of D matrix can be taken depending on the plane stress or the plane strain condition. In this analysis plane 
stress condition is considered for the deformation of the mild steel plate under loading. 
Now , by the finite element analysis method, 

KQ = F (5) 
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where, K is the global stiffness matrix, Q is the displacement matrix, and F is the force vector. 

Equation (3) is solved for the vector matrix Q which will give the values of displacement under the influence of the force 
vector F. 

III. PROBLEM STATEMENT 

FIG. 3 clearly elaborates the direction and the magnitude of the point loads acting on the plate . The plate is 20 mm in height 
and 30 mm in width and 10 mm in thickness. A point load of 100 N is acting in the negative Y direction at the free end of the 
plate and point load of 50 N in the positive X direction. The plate is fixed along the edge AB. 
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FIG. 3 

Now this triangular plate is mapped to a isoparametric constant strain triangle. The problem is solved using FEA 
and one element. Thus the points 1, 2 and 3 represent the nodes of the constant strain triangular element. Coordinate system 
used indicating the positive directions is shown in the figure. Point B is considered as origin for the simplicity of analysis. 

3.1) Material Properties of mild steel: 

The material properties of mild steel are shown in the following TABLE 1. 



TABLE 1 



Material 
Used 


Young's 
modulus 

(E) 


Poisson's 
Ratio (v) 


Density 
(P) 


Mild Steel 


210 GPa 


0.3 


7850 
kg/m 3 



IV. CALCULATIONS 

To study the deformation in the plate, under given force vectors and the dimensions of the plate , plane stress 
condition is assumed. Thus the value of matrix D is given as 

f \ f \ 
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Using eq. 4; 

det J = 2 x A = 2 x (- x20x30) = 600 
2 

Using eq. 3 , value of matrix B is found out. 
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k e = t e A e B T DB = 10 x 300 x B T DB 

Using the program MATLAB above complex multiplication has been carried out and the value for the element stiffness 
matrix is 
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In this case as shown in FIG 3; only one CST element is considered for analysis. Thus, 
k e =k g , where k g is the global stiffness matrix. 
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Using equation (5); 



Global force vector will be given bv s 
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.(6) 



The edge AB of the mild steel plate is fixed as shown in FIG .3 Thus the DOFs q l , q 2 , q 5 and q 6 will be 

constrained to be zero. So now using the elimination method to solve the equation no. (6) striking off the rows and columns 
corresponding to the constrained Degrees of freedoms, we get, 
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Solving the above equation , we get, 
q 3 = 6.5 x 10~ 5 mm ; 



q 4 =-3.71x 10 4 mm 



(7) 



Let 0" be the stress generated in the element. 

a = DBq , where q is the element nodal displacement vector. 
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a = {0.49, 0. 15, -0.99} r MP a 



(8) 
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V. ANSYS Modeling 

Commercially available software ANSYS 13.0 is used to model the problem and solve for deformation. ANSYS 
parametric design language (APDL) is utilized for this purpose. A model is generated in the work plane for the given 
dimension and thickness of the mild steel plate. After specifying the material properties and real constants as shown in 
TABLE 1, boundary condition is applied. That is both the nodes 1 and 3 are fixed thus constraining their displacement to 
zero. At the node 2, constant point load of 50 N in the positive X axis and 100 N in the negative Y axis direction is applied. 
The generated model is shown in FIG. 4. Deformation analysis finding the nodal displacements using the numeric method is 
done by the ANSYS postprocessor. The stress generated in the element is also found out. The displacement in the positive X 
direction of the node 2 is 0.65E-04 mm and in the negative Y direction 0.37E-03 mm. Using the command plot results- a 
deformed shape of the plate is generated by the ANSYS software. This is shown in the FIG 5 



FIG. 4 FIG 5 

VI. Results & Discussions 

From eq. (7) and (8) and the results obtained from ANSYS modeling, it can be conclusively averred that the results 
of the finite element analysis using Constant Strain Triangle and that of ANSYS show good agreement. The deformation of 
the plate subjected to given boundary condition and loading and the corresponding Von-Mises Stress generated were equal 
in both the methods. 

VII. Conclusions 

This paper mainly focused on finding the deformation of the triangular mild steel plate subjected to a boundary 
condition under the influence of point loads. The deformation analysis is first carried out by finite element method using the 
three node Constant Strain Triangle (CST) as a basic geometric shape. Later, for the same structure and load, boundary 
conditions, analysis has been performed using analysis software ANSYS (APDL). Finally, the results obtained from FEA 
and ANSYS are compared and they are closely converging. Thus use of Constant Strain Triangle as a master element for 

FEA in complex stress and deformation analysis is proposed because of the following reasons: 

1) The stress value in the basic modeling element is constant and depends on nodal displacements. This saves a lot of 
computational effort and time required for analysis. The shape functions are linear over the element. Isoparametric 
representation allows use of same shape functions for representing geometry and displacements inside the element. This 
approach lends simplicity of development and retains uniformity with other complex elements. 

2) The results obtained by FEA using CST as modeling element are in good agreement with the solution. 
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